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Abstract 

We simulate 2-flavour lattice QCD at finite isospin chemical potential fij, for temperatures 
. close to the finite temperature transition from hadronic matter to a quark-gluon plasma. The 

^ , fii dependence of the transition coupling is observed and used to estimate the decrease in the 

O ■ . o 

' transition temperature with increasing fij. These simulations are performed on an 8 x 4 lattice 

o ■ 

■rj" ■ at 3 different quark masses. Our estimate of the magnitude of the fluctuations of the phase of the 



o 



X 



fermion determinant at small quark-number chemical potential fx, suggest that the position of the 
small and small fij transitions should be the same for /i/ = 2/1, and we argue that the nature of 
these transitions should be the same. For all m < m.,^ the smoothness of these transitions and the 
values of the Binder cumulant B4, indicate that these transitions are mere crossovers, and show no 
sign of the expected critical endpoint. For fij > and a small isospin {I3) breaking term A, we 
do find evidence of a critical endpoint which would indicate that, for A = 0, there is a tricritical 
point on the phase boundary where the pion condensate evaporates, where this phase transition 
changes from second to first order. 
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I. INTRODUCTION 



QCD at finite baryon-/quark-number density describes nuclear matter. Beyond nuclei it 
describes the physics of neutron stars and has the potential to predict such exotic objects 
as quark stars. Hot hadronic matter at low baryon-number density was present in the early 
universe. Relativistic heavy-ion collisions at RHIC and CERN produce hot nuclear matter. 

QCD at a finite chemical potential fi for quark-number, has a complex fermion deter- 
minant, which makes the naive application of standard lattice simulation methods, which 
are based on importance sampling, difficult if not impossible. To circumvent these prob- 
lems people have introduced various schemes which are at 
small a. These include various reweierhtinsr techniques IfJ, 

phy*a. ob.e„ab.es a. powe. se.e. JflBSaflB or related pa.a„,ete.s M- 



icable to high temperatures and 
1^, and methods which exp and 



Another approach is to study theories which are expected to possess some of the properties 
of QCD at finite fi, but have real positive fermion determinants, making them amenable to 
standard simulation methods. One such theory is QCD at finite chemical potential ni for 
isospin (/a). This theory has been studied both by effective (chiral) Lagrangian techniques 
12, [l^, as well as by direct lattice simulations jl4|. At zero temperature these studies 



indicate that this theory undergoes a second-order phase transition with mean-field critical 
exponents at fij = rriT^, to a state characterized by a charged pion condensate which breaks 
I3 spontaneously. 

We report here a study of 2-fiavour lattice QCD at finite fii and finite temperature (T), 
in the neighbourhood of the finite temperature transition from hadronic matter to a quark- 
gluon plasma. Since we work at finite quark mass to make the pion massive and thus to move 
pion condensation to finite fij, the finite temperature transitions form a line of crossovers em- 
anating from the fij = transition, for small /i/. We calculate the position of this crossover 
as a function of fij on an 8^ x 4 lattice for 3 different quark masses (m = 0.05,0.1,0.2), 
from the peaks of the susceptibilities of the various observables, using Ferrenberg-Swendsen 
reweighting to interpolate between the /? values used in our simulations. For /i/ < m^, we set 
the symmetry breaking parameter A = 0. From estimations of the fluctuations of the phase 
of the fermion determinant for small quark-number chemical potential fi we shall argue that 
there is an appreciable range of fi over which these fluctuations are small enough that the 
position of the crossover at finite fi will be the same as that at finite fij with fij = 2fi. We 
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find good agreement with the /i dependence of this transition observed by de Forcrand and 
Phihpsen j6[. This agreement between the jj, and yU/ dependence of the transition (3 = 
and hence temperature was noted by the Bielefeld-Swansea group ^ . We also find that the 
transition for each of our 3 masses appears to remain a crossover with no sign of a critical 
endpoint for al l p^j < ttk- Preliminary results from these simulations have been presented 
at conferences UQ, 

We have also studied the finite temperature transition for yU/ > m^. Here, for symme- 
try breaking parameter A = 0, the pion condensate evaporates at the finite temperature 
transition, which is thus a true phase transition. However, since here the fermion propa- 
gator becomes singular (at least in the infinite volume limit) for temperatures below this 
transition, because of the Goldstone boson associated with the spontaneous breaking of J3, 
we are forced to work at finite (small) A, where the transition again becomes a crossover. 
Here we shall present evidence for a critical endpoint beyond which the transition becomes 
first order. Because A is small, we shall argue that this first order behaviour persists to 
A = 0. At A = 0, the finite temperature crossover is replaced by a second order transition, 
the first order transition remains first order and the critical endpoint becomes a tricritical 
point. Although most of our simulations were performed on 8^ x 4 lattices, we performed 
some simulations on 16^ x 4 lattices close to the critical endpoint. 

In section 2, we introduce lattice QCD at finite /i/. In section 3 we define the fourth order 
Binder cumulants which we use to study the nature of the transitions. Section 4 describes 
our simulations and results for small ni (/// < m^). The large ni simulations and results 
are presented in section 5. Section 6 contains discussions and conclusions. 



II. LATTICE QCD AT FINITE ^u/ 

The staggered quark action for lattice QCD at finite chemical potential /i/ for isospin 
(^3) is 



sites 



(1) 



where Jp^^r^fij) is the standard staggered quark transcription of ^with the links in the +t 
direction multiplied by exp(|r3/i/) and those in the —t direction multiplied by exp(— |r3/i/). 
The term proportional to A is an explicit Is = symmetry breaking term. This term serves 
two purposes. Firstly, such a term is necessary if one is to see evidence for spontaneous /a 



breaking on a finite lattice. Secondly, it prevents the Dirac operator from becoming singular, 
as we see below, ri, T2 and T3 are the isospin matrices so that this Dirac operator is a 2 x 2 
matrix in isospin space. The determinant 

detiK^rg/xj) + m + zAers] = det[^U + A^], (2) 

where 

A = p{^fii)+m, (3) 

is a 1 X 1 matrix in isospin space, which means that we only need use a single flavour- 
component fermion field in our simulations. This determinant is real and positive allowing 
us to use standard hybrid molecular-dynamics simulations, with noisy fermions to allow us 
to tune the number of flavours from 8 down to 2. 

We note that, for A = 0, the determinant of equation |2l is just the magnitude of the 
determinant for 8-flavour lattice QCD with quark-number chemical potential 

= ^/^/- (4) 
Observables for this theory include the chiral condensate, 

i'iPij) ^ ixx), (5) 

the charged pion condensate 

i{'4^1^r2ip) <^ i{xeT2x) (6) 

and the isospin density 

We will also be interested in the Wilson Line (Polyakov Loop), and the plaquette observable 

PLAQUETTE = 5n = 1 - -ReTr?7n. (8) 

3 

III. FOURTH-ORDER BINDER CUMULANTS 

If X is an observable, its 4-th order Binder cumulant is defined by 

_ ((X - {x)Y) 
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which approaches a universal vakie at a critical point It has been pointed out if one 
chooses X to be an eigenvector of the critical scaling Hamiltonian, i?4 will be as close as 
is possible to its infinite volume limit on finite volumes jl^. If one plots as a function 
of those parameters which parameterize the departure from the critical point, the curves 
obtained for different lattice sizes will intersect at the critical point if X is indeed an eigen- 
vector. For other choices of X the intersections of such curves will only tend to this unique 
value in the infinite volume limit. The value of the cumulant at this point of intersection 
will be that characteristic of the universality class of this critical point and the nature of 
the observable. 

For transitions other than critical points, the Binder cumulant only attains its charac- 
teristic value in the infinite volume limit. For a crossover, the infinite volume value for 
the Binder cumulant for the order parameter is -B4 = 3. For a first order transition, this 
Binder cumulant is i?4 = 1. The critical endpoint we are seeking is expected to be in 
the universality class of the 3-dimensional Ising model for which = 1.604(1). For a 
mean field critical point B4 = r(5/4)r(l/4)/r(3/4)^ = 2.1884... for a 1-component or- 
der parameter Q], or i?4 = 7r/2 = 1.570796... for a 2 component order parameter. At 
a 3-dimensional tricritical point for a 1-component order parameter i?4 = 2 [2^, and 
-B4 = r(l/3)/r(2/3)^ = 1.460998... for a 2-component order parameter. 



IV. SIMULATIONS AND RESULTS FOR m < m„ 

We have simulated 2-flavour QCD on 8'^ x 4 lattices in the neighbourhood of the finite 
temperature transition from hadronic matter to a quark-gluon plasma, for small values of 
the isospin (J3) chemical potential fij. Here small fii means fij < so that, even at 
zero temperature, the system is in the normal phase, i.e. in the phase where there is no 
Js-breaking charged pion condensate. We set A = 0, since a finite A is only needed when 
there is a possibility of spontaneous breaking. Simulations were performed at 3 different 
quark masses m = 0.05, 0.1, 0.2. 

At the lowest quark mass m = 0.05 we performed simulations over a range of fij values 
^ fJ'i ^ 0.55, where the highest fij value is only just below the critical fij above which a 
pion condensate forms at low temperature. The larger quark masses were chosen to allow 
an even larger range of /i/s for the normal phase. For m = 0.1 we simulated over the range 
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< /U/ < 0.7, while for m = 0.2 we covered the range < /x/ < 1. At each of the selected 
fii values we performed simulations over a range of (3 = Q/g"^ values spanning the finite 
temperature transition. We used a range of updating increments (in molecular-dynamics 
'time') for these simulations. These ranged from dt = 0.1 for m = 0.1,0.2 and small fij 
down to dt = 0.01 for m = 0.05 and /i/ = 0.55 close to the transition. For each quark mass 
we performed simulations as long as 20, 000 molecular-dynamics time units (trajectories) at 
individual values of (/3, fij) close to a transition. 

At each value of m, /i/ and P we measured the average plaquette, the Wilson Line 
(Polyakov Loop), the chiral condensate and the isospin density for each trajectory. For the 
fermionic quantities, where we calculate stochastic estimates, we used 5 noise vectors for 
each trajectory, which enabled us to make unbiased estimates of the susceptibilities and 
Binder cumulants. Figure ^ shows the Wilson Line as a function of (3 for each /i/ at quark 
mass m = 0.05. Note that there is a rapid crossover marking the transition. In addition we 
see that the position of the crossover moves towards smaller f3 and hence lower temperature 
as /i/ is increased. However, we notice that the crossover (3, (3c varies only slowly with 
/i/. The corresponding values of the chiral condensate, {ipip) are given in figure El Again 
we see a rapid crossover close to the position of that for the Wilson Line. Figure El shows 
the behaviour of the isospin density j'q for the same mass and /i/s. Here we see the finite 
temperature transition again. We note that the value of the isospin density in the quark- 
gluon plasma (high /?) increases with increasing /x/. At each /x/ it appears to level off at large 
j3. Note that this is not the lattice artifact of saturation; j'q in this domain is far below its 
saturation value of 3. The rise in jg occurs because increasing /x/ raises the Fermi surface. 
These observables for m = 0.1 and m = 0.2 behave very similarly to those for m = 0.05, 
except that the crossovers occur at larger (3 values as mass is increased. 

The transitions we have observed in each of these masses and chemical potentials appear 
to be smooth crossovers rather than actual phase transitions, as is expected to be the case 
for fij = (we will present further evidence for this later in this section). The position 
of the transition is thus defined as the (3 value which maximizes a chosen susceptibility. 
(Such definitions and Ferrenberg-Swendsen reweighting are used by other groups, as are the 
Binder cumulant methods used below P. Isl. IgI Rligjl.lThis is a reasonable definition only if 
the positions of the maxima of the susceptibilities for the various observables are close, at 
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SU(3) Nf=2 m=0.05 A=0 8^x4 lattice 
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FIG. 1: Wilson Line as a function of (3 for various /x/ < and m = 0.05. 



least in the infinite volume limit. The susceptibility for a chosen observable O is defined as 



Xo = V{0^ - {Of), 



(10) 



Where V is the space-time volume of the lattice. Note that this is the correct definition 
only for a local observable O. We have also used this definition for the Wilson Loop which 
is only local in 3-space. Thus what we call Xwuson is strictly NtXwuson = Xwuson/T. 

Figure m shows the plaquette susceptibilities for m = 0.05 at those /3s at which we per- 
formed our simulations for each /i/. These susceptibilities are clearly strongly peaked, and 
the peaks move to lower /3s as fij is increased. Figure El gives the corresponding suscepti- 
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SU(3) Nf=2 m=0.05 X=0 8^x4 lattice 



1.25 



1.00 - 



0.75 



0.50 - 



0.25 - 



0.00 



X 

o 
□ 
o 



- + 



Ml 
Ml 
Ml 

Ml 
Mi 
Mi 




J I I L 



J I I L 



J I L 



5.0 



5.2 



5.4 



FIG. 2: {ipip) as a function of (5 for various ^/ < and m = 0.05. 

bilities for the Wilson/Polyakov Line. Again these susceptibihties are strongly peaked and 
the peak moves to lower (3 as /x/ is increased. The main difference is that the peaks in these 
susceptibilities decrease slightly in height as yU/ is increased, whereas those for the plaquette 
susceptibilities increase with increasing /i/. Figure El shows the susceptibilities for the chiral 
condensate, also for m = 0.05. By using all 5 stochastic estimators of (ipip) and removing 
the noise-diagonal contribution we obtain an unbiased estimate of this susceptibility These 
susceptibilities are the most strongly peaked of all the susceptibilities and the height of 
the peaks increases with increasing hj. Finally we show the susceptibilities for the isospin 
densities at m = 0.05 in figure Again we used the 5 noise vectors to obtain an unbiased 
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SU(3) Nf=2 m=0.05 X=0 8^x4 lattice 
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FIG. 3: Jo as a function of /3 for various /i/ < and m = 0.05. 



estimator. The rapid increase in the height of these peaks with increasing fij corresponds 
to the increase in Jq seen in figure El 

In order to pinpoint the susceptibility peaks more precisely, we use the distribution of 
observables and plaquette actions measured during our runs and use Ferrenberg-Swendsen 
reweighting |2l| to estimate the susceptibilities at (3 values close to those at which we have 
performed simulations. If O is an observable for which Oi, i = l,...n are the measured 
values (lattice averaged), and Sni are the corresponding plaquette actions, at /? = Pq, then 

Erexp{-QV{(3-(3o)Sm)0, 



Eiexp(-6y(/5-/3o)5n,) 



(11) 



9 



SU(3) Nf=2 m=0.05 X=0 8^x4 lattice 
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FIG. 4: Plaquette susceptibilities as functions of /3 for m < m^, for m = 0.05. 

for P close enough to /?o that the distributions of Sn values at /? and /?o have significant over- 
lap. Applying this formula to estimate both (O) and (O^) yields the desired susceptibility 
Xo- Jackknife methods are used to determine the errors in both the susceptibilities and the 
positions of their peaks. It turns out that for our simulations at each m and fij, we have 
3 — 5 /? values close enough to the peak to be used to determine Pc- After checking that 
the estimates of Pc from each of these points are consistent, we obtain our final estimate as 
a weighted average of these. 

In figure |H1 we plot our Pc estimates from each of the 4 susceptibilities as functions of 
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SU(3) Nf=2 m=0.05 X=0 8^x4 lattice 
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FIG. 5: Wilson line susceptibilities as functions of /? for m < m,r) for m = 0.05. 
fi'j for all 3 masses. There is clearly good agreement between the (3c values obtained from 



the different suscepti 
references l22L ]2 



jilities. Arguments as to why this should be so have been presented in 
Since for a fixed mass 



= Pc{-l^i) (12) 

is a function of fij. For \fii\ small enough it is also an analytic function of fij. Thus for 
small /i/, 

Pcifii) = a + b^j + cfij + ... (13) 
We have therefore fit Pc to the form Pdf^i) = a + bfij, for each mass. Choosing to fit the 
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SU(3) Nf=2 m=0.05 X=0 8^x4 lattice 
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FIG. 6: ip^jj susceptibilities as functions of (3 for fij < m.,^, for m = 0.05. 

plaquette susceptibilities, we get 

/?c = 5.3198(9) - 0.134(6)/i^ form = 0.05 
pc = 5.3731(11) - 0.126(5)/i5 form = 0.1 
(3c = 5.4443(27) - 0.118(4)/i5 for m = 0.2. 



(14) 



These lines are plotted in figure |H1 Since x^/dof is 3.2 for m = 0.05, 1.1 for m = 0.1 and 
2.4 for m = 0.2, these are not excellent fits. However, including the quartic term, does not 



12 
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FIG. 7: Isospin density susceptibilities as functions of /? for < m-,^, for m = 0.05. 



improve the m = 0.05 fit, and adding even more terms would be a meaningless exercise. 

It was observed by the Bielefeld-Swansea group that at small fi and fij, the dependence 
of on /i and on fij was identical within their errors |^]. In our notation this would mean 
that 

=/5c(/i/ = 2/i). (15) 

This means that one can determine the position of the finite temperature transition for small 
/i by simulating with the magnitude of the fermion determinant and ignoring the phase. Let 
us give an intuitive argument why this might be the case. If 6 is the phase of the fermion 
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FIG. 8: as functions of fi'j together with straight hne fits for each mass. The bottom set of 
points and hne are for m = 0.05. The middle set of points and hne are for m = 0.1. The top set 
of points and hne are for m = 0.2. 

determinant, and (9 is a gauge- field observable (such as the plaquette or Wilson line) then 

(exp(z^)0)^,=2M 



(16) 



(exp(i^))^,=2/. 

where for O real we can replace exp{id) with cos(6'). For the lattice size used by the 
Bielefeld-Swansea group, there was a range of fi over which this denominator was not too 
small and varied smoothly with /i and slowly with (3, as f3 swept across the f3c- Here the 
denominator cannot be responsible for the transition. When the fluctuations in 6 are so well 
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behaved, it is reasonable to treat exp{i6)0 as another observable. Since the position of the 
transition appears to be nearly independent of the chosen observable, this would suggest 
that the position of the transition would be the same for this new observable. If so, the 
smooth behaviour of the denominator would imply that the position of the transition at 
finite fi should be the same as that for the transition at finite fij for fij = 2/i. It is also 
not unreasonable to assume that the nature of the 2 transitions might be the same. Such 
observations are not new (see for example [2^). 

If the relation between /?c(/^) and PdfJ'i) holds with the standard staggered action (the 
Bielefeld-Swansea group used the p-4 action) we can compare our formulae for Pdf^i) (equa- 
tion [T3|) with that obtained by de Forcrand and Philipsen 

Pc = 5.2865(18) - 0.149(10)/i? for m = 0.025 (17) 

where we have made the substitution fi = /i//2 in their equation. This would appear to 
be consistent with our equations, taking into account the difference in mass. To examine 
whether this agreement is quantitative, we fit equations ll4fl7l to the expected scaling form 

Pc{m,fij) = ^c(m, 0) + a(m)/Xj 
f3c{m,0) = /?c(0, 0) + 6m/5^ 

a(m) = a{0)+CTn^. (18) 



Such scaling fits have been considered by 27|, |28[ at zero chemical potentials. For the 



expected continuum 0(4) scaling l/f3mS ~ 0.55, while for the lattice 0(2) scaling l/(3mS ~ 
0.59. [Note that such scaling is only derivable for the case of finite /i, where, in the chiral 
limit, the line of crossovers becomes a line of second order transitions in the same universality 
class as the fi = 0(/i/ = 0) transition. We are using the assumed relationship between finite 
fi and finite fij to extend it to finite fij.] The fit of all 4 equations to 0(4) scaling gives 
/3c(0,0) = 5.210(3), b = 0.57(1), with a x^/dof = 1.6, and a(0) = 0.152(6), c = 0.85(19) 
with a x^/dof = 0.5. The fit to 0(2) scaling gives pc{0,0) = 5.219(3), b = 0.59(2) with 
a x^/dof = 2.4 and a(0) = 0.151(6), c = 0.087(20) with a x^/dof = 0.5. Considering 
the quality of the fits in equations El we consider either of these scaling fits to be good 
enough to support our claim that we are consistent with de Forcrand and Philipsen, and 
that the combined measurements are consistent with the expected scaling with quark mass 



m. Note that our value of /5c(0, 0) is less than that obtained in |27|, |28|. We should not 
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expect good agreement with the later paper, since it uses a larger lattice and finite size 
effects are non-negligible on an 8^ x 4 lattice. The fit in the earlier work was over the mass 
range 0.02 < m < 0.075, while ours was over the range 0.025 < m < 0.2. Considering the 
rapid variation of the scaling form at small m, the difference between our result, 5.210(3) 
and their's, 5.222(3) is perhaps not surprising. 

What remains to be checked is that the phase (6) of the fermion determinant is well 
behaved. Since calculating the fermion determinant is very expensive, we use the series 
expansion for 9 given in In our normalization, 

9 = ^/i/yim(jo) + 0(/i?) (19) 

where jo is the number density normalized to 4 flavours (1 staggered fermion fleld). We 
use our 5 stochastic estimator s/conflguration of jo to obtain an unbiased estimator of {9"^) 
through order fij. (We also made an unbiased estimator of (jo) which had a poor enough 
signal/noise ratio that we did not even try to estimate (0^) or the higher order contributions 
to {9"^)-) Our results for a range of j3 values which span the fij = transition for each 
quark mass are given in table IH A reasonable measure of how "well-behaved" this phase 
is, is (cos(^)). When this quantity is close to 1, the oscillations in phase are small, and it 
is reasonable to produce ensembles with the magnitude of the determinant and to include 
the phase in the measurement. When this expectation value falls towards zero 9 is almost 
uniformly distributed over the interval (— tt, +7r], and the contributions of conflgurations 
generated using the magnitude of the determinant can easily cancel, as they would in the 
denominator of equation for this case. How small {cos{9)) can get before generating 
ensembles without the phase becomes invalid is a matter of "experimentation", but one 
might expect that (cos(^^)) > 0.5 would be a reasonable range over which we could use this 
method. To the order in fij to which we work we must take cos(6') ^ 1 — ^9^. Applying 
our criterion to the measurements of table U we see that, even in the worst case we should 
be able to trust the relationship between measurements at flnite fij and flnite /i, out to 
/Xj ~ 0.15, i.e. out to fij ~ 0.4. Thus it is not unreasonable to assume that this relation will 
be a reasonable approximation for most of the region /i/ < m,r- 

Let us now examine the nature of these transitions more closely. We have observed that 
the transitions appear smooth in all measured observables. This suggests that they are 
merely rapid crossovers. Histogramming those observables which could show discontinuities 
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TABLE I: Fluctuations in the phase of the fermion determinant. 

if there were a first order transition, shows a single broad peak for all masses considered for 
all fij < rriTj, which suggests a crossover (or possibly a second order transition) but not a 
first order transition. We note that on such small lattices, one can observe a double peak 
structure, even where the transition is a crossover or second order transition. However, it 
is rare that a first order transition would not show a double peak, unless it were very weak. 
In figure El we show histograms of the Wilson Line (Polyakov Loop) for m = 0.05 at an 
intermediate value of fij (0.3) and one close to (0.55). These both show a single broad 
peak as advertised, and are typical. [We chose to show the Wilson Line rather than the chiral 
condensate, since use of stochastic estimators (even after averaging over all 5 estimates for 
each configuration) could possibly obscure a double peak.] 

Finally, we have calculated the fourth order Binder cumulants B4 for the chiral condensate 



17 



8^x4 LATTICE iS = 5.30?5 /ii^ = 0.3 
600 I ^ 1 




0.0 0.2 0.4 0.6 

WILSON LINE 



8^x4 LATTICE |8=5.2825 /Xi = 0.55 
500 I 1 




WILSON LINE 

FIG. 9: Histograms of distribution of Wilson Line values for m = 0.05: a) For /i/ = 0.3, (3 = 5.3075; 
b) For = 0.55, /? = 5.2825. 
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at the transition for each m and Having 5 noisy estimators per configuration, we were able 
to generate an unbiased estimator for B4. We again use Ferrenberg-Swendsen reweighting 
to interpolate between those (3 values at which we ran our simulations. We determined the 
position of the transition for each m and fij as that j3 which minimized -B4. This method of 
determining the position of the transition gave f3c values in excellent agreement with those 
obtained from the maxima of the corresponding susceptibilities. We plot these values in 
figure Uni 
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FIG. 10: Fourth order Binder cumulants {B4) for (ipip) as a function of fij. The dashed line is at 
B4 = 1.604, the value for the 3-dimensional Ising model. 
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If there were a critical endpoint, which is expected to be in the Ising universahty class, 
we would expect -B4 to pass through its Ising value, B4 = 1.604, at this endpoint. For 
the crossover region, B4 should lie above the Ising value, approaching 3 in the limit of 
large lattices. In the first order domain (if it existed) B4 should lie below the Ising value, 
approaching 1 in the large lattice limit. We have plotted the Ising value as a dashed line 
in figure [TUl Clearly B4 lies well above 1.604 and shows no sign of approaching this value. 
Hence the evidence from Binder cumulants supports our conclusion that the transition is a 
crossover for all fij < m^, over a range of quark masses. 

V. SIMULATIONS AND RESULTS FOR /x/ > 

In the region where fij > m-Tr, for A = 0, the charged pion condensate evaporates at 
the finite temperature transition and I3 symmetry, which is broken spontaneously at low 
temperature, is restored. Hence this finite temperature transition is a true phase transition. 
However, since in this case the low temperature phase has a true Goldstone mode, this would 
render the Dirac operator singular (at least in the large volume limit). Hence we must use 
a non-zero A, which we keep small so that we can infer information about the A — limit. 
For A 7^ 0, the phase transition is no longer required, and our earlier work indicates that for 
Hi just above m.,^, the transition is softened to a crossover. We can now search for a critical 
endpoint in the high fij {fii > uIt,) regime. Again the critical endpoint would be expected 
to lie in the universality class of the 3-dimensional Ising model. For /i/ above this endpoint 
the transition would be first order. Unfortunately, in this domain, we cannot argue that the 
finite /// and the finite /i transitions are related. 

We have performed simulations on an 8^x4 lattice at quark mass m = 0.05 and A = 0.005, 
at /i/ = 0.6, fij = 0.8 and /i/ = 1.0. The /i/ = 0.8 simulations were repeated on a 16^ x 4 
lattice. In figure ^2 we show the behaviour of the Wilson Lines as functions of f3 for the 3 
fij values from our simulations on an 8^ x 4 lattice. At fij = 1.0 we have only performed 
simulations very close to the transition. We obtained the high statistics needed to reveal the 
true nature of this transition at those (3 values closest to the transition for each /i/ — for 
/i/ = 0.6 we obtained 40,000 time units aX (3 = 5.27, for fij = 0.8 we obtained 40,000 time 
units at /5 = 5.265 using dt = 0.02 and a further 40,000 time units using dt = 0.01, while at 
/i/ = 1.0 we obtained 40,000 time units at /5 = 5.263. We show histograms of the Wilson 
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FIG. 11: Wilson Line (Polyakov Loop) for large /i/ on an 8'^ x 4 lattice with m = 0.05 and A = 0.005 

Line at these (3 and /i/ values in figure IT^ At yU/ = 0.6, the histogram shows no structure 
to suggest anything but a crossover, which would then become a second-order transition as 
A ^ 0. By /i/ = 0.8 we begin to see clear signs of double peak suggestive of a 2-state signal. 
The signs of a double peak and a 2-state signal persist at yU/ = 1.0. 

Since lattices as small as 8^ x 4 can show signs of a 2-state signal at a second order 
transition or even a crossover, we need to examine this transition more closely. For this 
reason we have performed simulations on a 16^ x 4 lattice at /i/ = 0.8. Figure ^1 shows 
the Wilson Line and the pion condensate from these simulations as functions of (3. The 
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FIG. 12: Wilson Line histograms for large /x/ on an 8^ X 4 lattice with m = 0.05 and A = 0.005. 
a) nj = 0.6, p = 5.27; b) m = 0.8, /3 = 5.265, dt = 0.1; c) m = 1.0, /3 = 5.263. 

reason dt was decreased from 0.02 to 0.01 close to the transition was that finite dt effects 
at dt = 0.02, both here and in our 8'^ x 4 runs at the same iij, can artificially enhance the 
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FIG. 13: a) Wilson Line as a function of /? at /i/ = 0.8 on a 16^ x 4 lattice, b) Charged pion 
condensate as a function of /3 at /// = 0.8 on a 16^ x 4 lattice. 
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2-state signal, dt = 0.01 appears free from such enhancements. For our dt = 0.01 runs at 
P = 5.263, P = 5.264 and /5 = 5.265 we ran for 30,000 time-units per /3 to obtain adequate 
statistics (for /? = 5.266 and /? = 5.267 we ran for 20,000 time units per /?). Figure IT^ shows 
a histogram of the Wilson Line values from our (3 = 5.265 runs. Although there is some 
double peak structure, the peaks are considerably closer together than they were for the 
8'^ X 4 lattice, a sign that the double peak structure is a finite volume artifact and not the 
sign of a true 2-state signal indicating a first order transition. 

16^x4 LATTICE )8=5.265 yUi=0.8 X=0.005 
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FIG. 14: Histogram of Wilson line values close to the transition on a 16^ x 4 lattice at fij = 0.8 
(/3 = 5.265). 

To clarify the situation we again turn to fourth order Binder cumulants. Here the obvious 
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choice is to look at the Binder cumulants of the pion condensate, since this is the order 
parameter of this transition in the A — > hmit. We plot versus fij in figure IT31 obtained 
using Ferrenberg-Swendsen reweighting to obtain B4 at that f3 which minimizes B4 for that 
particular value of fij. The 8^ x 4 points suggest that the Binder cumulant crosses the Ising 
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FIG. 15: Binder cumulants of the charged pion condensate, i{'4''y5T2ip) , as functions of fii 



value somewhere above fij = 0.8 and probably close to fij = 1.0. If so, this would indicate 
that there is a critical end point with Ising critical exponents at /i/ = fic with fic ~ 1.0. For 
fii > fic the transition would become first order. The 16^ x 4 Binder cumulant at fij = 0.8 
is large enough to indicate that /ic is indeed greater than 0.8. We would expect that if the 
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transition is first order for A = 0.005 it will also be first order for A = 0. Hence there will 
be a tricritical point for yU/ = fit where < A^c- Since A = 0.005 is rather small, we expect 
that fj,t ~ fj,c- 



VI. CONCLUSIONS 

We have simulated lattice QCD with 2 flavours of staggered quarks ('half a staggered 
quark field) with a chemical potential /i/ for isospin {I3), in the neighbourhood of the finite 
temperature transition. For fij < we have determined the /i/ dependence of Pc, the 
transition (3 = Q/g"^, for each of 3 quark masses. We have noted that the fluctuations of the 
phase of the fermion determinant on an 8^ x 4 lattice are well enough behaved for small /i/ 
that there should be a range of fij for which the dependence of /5c, and hence temperature, 
on /i/ and on the quark-number chemical potential fi should be identical for hj = 2fi, as was 
observed previously by the Bielefeld- Swansea group for the p-4 action Qj. 

What we find is that jSc falls slowly with increasing fij. This falloff is approximately linear 
in fij over the whole fij < region. The value of (3c a.t = increases with mass and 
the falloff with increasing fi'j becomes less steep as the mass is increased. This dependence 
on mass is small. We have taken the results of de Forcrand and Philipsen j^] and converted 
them from a /i dependence to a fij dependence. Since these were calculated at a smaller 
mass m = 0.025 than ours, we cannot compare them directly. We fit the mass dependence 
of our 'data' and theirs (a total of 4 different masses) to that expected from critical scaling 
with both the continuum 0(4), and lattice 0(2) universality classes, as was done at zero 



chemical potential in [27|,|28|. Both these fits prove acceptable. This is a direct confirmation 
that the fii and fi dependence of the transition temperature are the same, at low chemical 
potentials. 

Since the phase of the fermion determinant is an extensive quantity, the fluctuations in 
this phase increase with volume. This would suggest that the relation between finite n and 
finite /i/ transitions would fail as the spatial volume is increased (the temporal extent is 
fixed at 1/T), which is disturbing since the infinite volume limit is of the most physical 
importance. We suggest that the relevant phase to consider in establishing this relationship 
is not that on an arbitrarily large lattice, but is rather the phase on a lattice whose spatial 
size is of the order of the correlation length. Then we could limit our considerations of phase 
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fluctuations to a modest lattice size, unless we were very close to a critical point. 

It is worthwhile quantifying what we mean when we say that the dependence on fij is slow, 
and what this means for the physical quantity, temperature. If we assume that the transition 
temperature at /ij = /x = is ~ 173 MeV, then for our m = 0.05 simulations, 2-loop 
running of the coupling would imply that by /i/ = 0.55 or in physical units fij ~ 362 MeV 
(yU ^ 181 MeV), the transition temperature will have fallen to Tc ~ 164 MeV. The relevance 
of this is even more clear when one considers that this latter /i value is an order of magnitude 
larger than those chemical potentials believed accessible by RHIC. 

The smoothness of the transitions for fij < m,r for all 3 masses, and the absence of any 
sign of a 2-state signal strongly suggests that there is no critical endpoint, beyond which 
the transition would become first order, in this domain. Analysis of the 4-th order Binder 
cumulant for each of the transitions, yields values ^ 2. Since this quantity should pass 
through the 3-d Ising value 1.604(1) at a critical endpoint and lie below this value in the 
first order region, this validates our assumption that the finite temperature transition from 
hadronic matter to a quark-gluon plasma remains a crossover throughout this region, and 
suggests that there is no critical endpoint for the corresponding range of fi. 

We have also studied fij > m.^, where the finite temperature transition for symmetry 
breaking parameter A = is a true phase transition from a pion condensed superfiuid to a 
quark-gluon plasma. Here we have performed simulations with a small A (0.1m), where the 
second order transition for fij just above mT^, softens to a crossover. We see evidence that for 
III sufficiently large (~ 1), there is a critical endpoint, where the 4-th order Binder cumulant 
passes through its 3-d Ising value, and beyond which it is first order. This observation 
needs to be confirmed on larger lattices, where the passage through the critical endpoint is 
expected be considerably more rapid. Unfortunately, in this regime, we cannot argue that 
a critical endpoint in fij is in any way related to a corresponding critical endpoint in fi. 

We have recently begun to extend this work to the 3-fiavour case. Not only is this more 
physical, but one can argue that it is possible to tune the critical endpoint to be as close to 
/i/ = as one desires, by careful choice of the quark mass m. In particular we can choose the 
critical value /i/ = fic to obey fic < ^tti and lie in the domain where the and /i/ transitions 
are related. Studies of the 3 and 2 + 1 fiavour transitions by various methods have located 
such a critical endpoint, but their predictions of its location are not in agreement 
Hence we have a chance to clarify the situation by a more direct approach. 
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It has been pointed out by de Forcrand, Kim and Takaishi 29|, that simulating with finite 
fii provides a better ensemble for reweighting methods for finite /i, than simulations with 
zero chemical potential. Combined with our observations, this suggests that such reweighting 
would be optimal close to the finite temperature transition for small fi, and could be expected 
to give good predictions for observables in this domain. This should enable us to determine 
the equation-of-state in this low-/i domain. Of course, such reweighting requires calculating 
the phase of the fermion determinant. Doing this precisely would be prohibitively expensive 
for all but the smallest lattices. New methods for approximating the fermion determinant 
show promise for making these reweighting methods practical (3fll |. 
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